MIXED MODELS WITH n > 1 AND LARGE 
SCALE STRUCTURE CONSTRAINTS 



S.A. Bonometto 

Dipartimento di Fisica deU'Universitd & INFN sez. di Milano 
9P ' Via Celoria 16, 

120133 Milano, ITALY 
e-mail: bonometto@mi.infn.it 



0^ 



^ ■ E. Pierpaoli 



Oh' 



O 



SISSA - International School for Advanced Studies, Via Beirut 2/4, 

134013 Trieste, ITALY 



> 

in . • ^ • 

^f) e-mail: pierpa& sissa.it 

o ■ 
o 

oo 



Abstract 



Recent data on CBR anisotropies show a Doppler peak higher than expected in 
CDM cosmological models, if the spectral index n = 1. However, CDM and LCDM 
^ \ models with n > 1 can hardly be consistent with LSS data. Mixed models, instead, 

' whose transfer function is naturally steeper because of free-streaming in the hot 

component, may become consistent with data if n > 1, when O/j is large. This is 
^ , confirmed by our detailed analysis, extended both to models with a hot component 

5h \ whose momentum space distribution had a thermal origin (like massive neutrinos), 

and to models with a non-cold component arising from heavier particle decay. In 
this work we systematically search models which fulfill all constraints which can be 
implemented at the linear level. We find that a stringent linear constraint arises from 
fitting the extra-power parameter T. Other significant constraints arise comparing 
the expected abundances of galaxy clusters and high-z systems with observational 
data. If low values of T are permitted, mixed models with 1 < n < 1.4 can have up 
to ~ 45% of non-cold component, without violating any further linear constraint. 
Keeping to models with F > 0.13, a suitable part of the space parameter still 
allows up to ~ 30% of hot component (it is worth outlining that our stringent 
criteria allow only models with 0.10^ 0.16, if n < 1). Spectra with n > 1 
are briefly discussed within the frame of inflationary theories. We also outline that 
models with such large non-cold component would ease the solution of the so-called 
haryon catastrophe in galaxy clusters. 
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1 Introduction 

We need three kinds of ingredients to define a cosmological model: the back- 
ground metric, the substance mix and the form of early departures from 
homogeneity. Background metric features are set by the density parameter 
Vt = pI Per and the Hubble parameter H = /i 100 km s^^Mpc^^ (p is the av- 
erage density and pcr = SH"^ /SttG); the cosmic substance is fixed by partial 
density parameters, like fib = Ph/ Per for baryons, fie = Pc/ Per for cold-dark- 
matter (CDM), or A for vacuum; finally, early deviations from homogeneity are 
described by the amplitude and the spectral index n of the initial fluctuation 
spectrum. 

However, while early spectral features essentially do not depend on background 
and substance parameters, the linear evolution of the spectrum does and this 

is usually expressed through transfer functions Tz{k) {k = 2n/L, L is the co- 
moving length-scale), which give the ratio between the spectrum at a redshift 
z and its initial value. 

In this paper we shall restrict ourselves to spatially flat models, with no resid- 
ual vacuum energy, and focus on the interplay between the nature of the 
cosmic substance (yielding %) and primeval spectral features, performing a 
systematic analysis of the amount and nature of hot dark matter (HDM), al- 
lowed by observational constraints, once we take the primeval spectral index 
(see below) in the range 1 < n < 1.4. 

In this work we shall however consider only large scale structure (LSS) fea- 
tures obtainable from linear theory, together with data on cosmic background 
radiation (CBR) anisotropics, which also refer to an epoch when fluctuations 

were still linear. The CBR temperature is almost direction independent and, 
once the dipolc term (which accounts for the motion of Earth in respect to 
CBR) is subtracted, its fluctuations AT{9,4>)/T = J2i,mC'imyim{(^,4>)^ 10~^. 
Since we expect that deviations from a Gaussian behaviour are substantially 
absent at these stages, anisotropy features are suitably expressed through the 
quantities Q = (|a;mp), for / > 2. As we shall sec, CBR and LSS linear con- 
straints already cause a severe selection among models and models to be used 
in numerical simulations should be submitted flrst to linear tests. 

Let us now consider recent data on CBR anisotropics. As usual, we shall 
assume that the primeval spectrum reads 

P^{k) ^ A^x^ixk)''-'^ (1.1) 
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{x = Xo-Xrec IS thc distaiicc from the recombination band, the present 

horizon radius). It is then known (see, e.g., Ma & Bertschinger 1996) that 

Ci = (4 Trf J dk k^P^{k) A^{k, Xo) (1.2) 

with 

Xq 

Ai{k,Xo) = j dxS{k,x)ji[{xo- x)k] (1.3) 



(Seljak & Zaldarriaga 1996). In this general expression, the source function 
S{k, x) contains contributions originating both at the last scattering band and 
along the path from it to us. For Q = 1 and low I (large scales L), the most 
relevant effect is Sachs- Wolfe, and the simple form 

Ci ~ {2n/3yA^ J {dk/k){xokr-'jf{x,k) (1.4) 

can be used, which does not depend on the substance characteristics, but only 
on n. Although this approximation holds for rather small / only, in most cases 
the nature of the substance does not drastically affect the / dependence of Ci 
up to the so-called Doppler peak (/ ~ 200) (Dodelson et al. 1996). 

Using eq. (1.4) and measuring Ci up to Z ~ 30, the COBE collaboration 
(Bennet et al. 1996) could estimate the spectral index and found n = 1.2 ± 
0.3. Therefore, COBE data already favour n > 1, although they are fully 
compatible with n = 1. 

More recent measurement (Netterfield et al., 1997), extending up to Z ~ 400 

(corresponding to angles ~ 0.5°), explored the doppler peak. The detected 
values seem to exclude models with low Q values (Hancock et al. 1998), and 
agree with CDM and fl = 1 only for h = 0.3 (Lineweaver 1997, Lineweaver & 
Barbosa 1998). An alternative possibility to fit such data is that n> 1. 

Although these results obviously need further confirm, they seem to tell us 
that inflation generated a spectrum P{k) with n > 1. 

There are two main reasons why this is seen as a problem by some cosmologists. 
In the first place, the generic prediction of inflation is that n ~ 1. Section 4 
will be devoted to a brief discussion on a class of inflationary theories which 
however predict n > 1. Furthermore, CDM fit with LSS improves for n < 1. 

Another way to improve the fit of models with LSS consists into assuming 
that DM has also a hot component. The basic philosophy of mixed dark 
matter (MDM) models, since the first computation of their transfer func- 
tion (Bonometto & Valdarnini, 1984; see also Achilli et al. 1985; Valdarnini & 
Bonometto 1985), is that the unexplored non-interacting world (gravitation 
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apart) can be quite complicated and probably contains several particle com- 
ponents. If one of them becomes non-relativistic when the scale L enters the 
horizon, then T is damped over all scales L. The detailed expression of T 
depends on the DM mix and the effect cannot be neglected if such component 
gives a substantial contribution to the overall density and L exceeds galaxy 
scales. 

An attempt was also performed to consider n > 1 in association to mixed DM 
(Lucchin et al. 1996), trying to balance opposite departures from standard 
CDM. However, the mix studied by these authors produces too many clusters, 
as was outlined by Borgani et al. (1996). Liddle et al. (1996) considered models 
with n > 1, in association with high Vth and possibly low h values, but they 
saw the same problems outlined by Borgani et al (1996) and did not indicate 
an area of fit with data. If it is confirmed that n > 1, however, there seems to 
be no alternative to try to follow such pattern. 

Owing to that, we tested a wide range of DM mixtures in association with 
values of n between 1 and 1.4. We found that the above difficulties can be 
overcame for a wide set of such models. Our analysis was restricted to the 
linear constraints described in sec. 3 and to models with Q = 1, Qf, = 0.1, 
A = and h — 0.5. Models were parametrized by Qh (the density parameter 
of the non-cold component) and d = z^q/ Zder'-, here we defined z^q = 10^; 
this would be the redshift when the photon density reaches the baryon+DM 
density, if all DM is CDM and 6 = r/2.733K = 1. The derelativization 
redshift Zder, instead, occurs when the average particle momentum (p) = m 
(particle mass). 

A particle model allowing a continuous selection of [r2/t,(i] pairs is obtained 
if heavier particles (X) decay into lighter massive particles {v), whose phase- 
space distribution turns out to be quite different from usual massive i/'s. Parti- 
cle models leading to such hot DM component were discussed in several papers 
{e.g., Bonometto, Gabbiani & Masiero 1994, Ghizzardi & Bonometto 1996) 
and the evolution of density fluctuations in such context was discussed by 
Pierpaoli & Bonometto (1995) and Pierpaoli et al. (1996). Hereafter we shall 
distinguish these models from models with phase-space particle distribution 
of thermal origin, calling volatile their non-cold component. However we shall 
show that also a large deal of models based on massive i/'s are consistent with 
linear constraints. 

Owing to the peculiar phase-space distribution of volatile DM and to the 
actual energy dependence of fluctuation amplitudes in such models, public 
algorithms cannot be used to evaluate T. Here we shall make use of a flexible 
algorithm that we developed also to this aim. A discussion of phase-space 
distributions, fluctuation evolution and constraints arising from particle kine- 
matics is given in sec. 2 and in Appendix A. More details on the algorithm are 
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given in Appendix B. As already mentioned, sections 3 and 4 deal with linear 
LSS constraints and inflationary theories yielding n > 1, respectively. Section 
5 will be devoted to the analysis of results concerning acceptable models. We 
provide some plots showing the allowed area in the parameter plane and some 
sample transferred spectra. Furthermore, for parameter choices consistent with 
hot DM {e.g. massive neutrinos), we compare the transferred spectra, finding 
small (expected) discrepancies between volatile models and those with massive 
i/'s. Section 6 is devoted to a final discussion and to conclusions. 



2 Pcirameter space limitation 

The substance of mixed models is made of various components, which can be 
classified according to their behaviour at the time when galactic mass-scales 
enter the horizon. Particles which become non-relativistic before then, are said 
cold. There can be one or more of such components, indistinguishable through 
astrophysical measures. Tepid or warm components were sometimes considered 
(see, e.g. Colombi et al. 1996; Pierpaoli et al., 1998), made of particles which 
become non-relativistic when small galactic scales (~ IO^Mq) have entered 
the horizon, but large galactic scales (~ IO^^Mq) still have not. Finally, hot 
component (s) are made of particles which become non-relativistic after the 
latter scale has entered the horizon. 

Neutrinos, if massive, are a typical hot component. They were dynamically 
coupled to radiation down to a temperature T^^dg ~ 900 keV. If their mass m <^ 
T^^dg, their number density, at any later time, is Ui, = {3({3)/4:iT^)g^T^ (after 
electron annihilation T,^ = T,^,dg[0'iy,dg / (^it)] = (4/ll)^/^T, where T is radiation 
temperature) and their momentum distribution (normalized to unity) reads 

3C(3)exp[p/r„(t)] + l '^"> 

also when p <^ m. Henceforth, when T <^ m, their distribution is not thermal, 
although its shape was originated in thermal equilibrium. Notice that, for high 
p, is cut off as exp(— p/Tj,). 

Using the distribution (2.1) we can evaluate 

(P) = (77rVl80C(3))T, = 3.1527; . (2.2) 

Accordingly, {p) — m. when T^, = 0.317 m and, for massive i/'s, 

d = 5.2972 • 4/iV(m/eV) , (2.3) 

while 

= 2.1437 • lO^V WeV)^V4^^ (2-4) 
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Hot components can also originate in the decay of a heavy particle X. Let 
Nx.dg be the comoving number density of X's at their decoupling, taking 
place at some early time t^g, well after they became non-relativistic, and let 
mx be their mass. At t ^ tdg their comoving number density reads: 

Nx{t) = Nx,dg exp[-(i - tdg) /ray] (2.5) 

with tdg <^ Tdy (decay time). We shall assume a two-body decay process, as is 
more likely for dynamical reasons: 

X^v + (p (2.6) 

The decay gives rise to a light {volatile) particle v, of mass m -C mx, and to 
a massless particle 0. It can be either a photon (7) or a sterile scalar, as is 
expected to exist in theories where a global invariance is broken below a suit- 
able energy scale (examples of such particles are familons and majorons) . We 
shall show that the latter case (sterile scalar) is potentially more significant. 

Once the decay process is over, the volatile momentum distribution reads: 

$„(p,i) = 2(Q/p)exp(-g) where Q = p'/p' (2.7) 

and 

P= {mx/2)[ady/a{t)] (2.7a) 

provided that X's, before they decay, never attain a density exceeding rela- 
tivistic components. As is shown in Appendix A, this is avoided if the condition 

nhd<^[l + 2/{l + g^/8.8)]-^ ior-f decay. (2.8) 

is fulfilled. Here refers to massless i/'s; e.g., for g^, = 6, eq. (2.14) yields 
Qhd < 0.46. If matter dominance occurs, the distribution (2.7) is distorted at 
small p, i.e. for particles originating in early decays. In this case the distribu- 
tion can still be evaluated numerically, but no simple analytical expression can 
be given. Particles with large p, instead, were born in decays occurring when 
radiation dominance was recovered and, therefore, the distribution is however 
cut off as exp(— p^/p^). The cut-off is therefore sharper than in the neutrino 
case. 

Using the distribution (2.7), it is easy to see that the average 

(P) = {V^/4:)mxady/a{t) = pV^/2 , (2.9) 

according to eq. (2.7a). v^s will therefore become non relativistic when p = 
2m/-^/7r; henceforth 

p= {2md/y^){z/Zeg) (2.10) 

(let us recall that we defined z^g = 10*^). Eq. (2.10) can be conveniently used 
in the distribution (2.7), instead of eq. (2.7a), as involves model parameters 
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(m and d), instead of early particle features {mx and Tdy), which should be 
further elaborated in order to work out aay — diTdy)- 

Let us now discuss the decay patterns in more detail. First of all, if = 7, 
the decay should occur at a redshift z'^ 10^, so to avoid distorsions of CBR 
spectrum (see, e.g. Burigana et al., 1991). Two further possibilities should 
then be considered. 

If the decay takes place before big-bang nucleosynthesis (BENS), the energy 
density of v''s is hmited by the requirement that BENS yields are consistent 
with observed nuclide abundances; henceforth the number of relativistic spin 
states 7 (see, e.g. Olive et al. 1997). Therefore, the actual limit on the 
energy density of v's depends on the number of relativistic i/'s at BENS. If 
some u has a mass exceeding ~ 1 MeV (particle data surely allow it for z/,-) 
and decays before decouphng, ^;'s can occupy its place during BENS. As far as 
cosmology is concerned, differences between such volatile DM and ordinary hot 
DM would arise from the following 3 kinematic reasons: (i) the different phase- 
space distribution of v's in respect to massive i/'s; (ii) the different energy 
dependence of fluctuation amplitudes; (iii) the allowed continuous range of 
HDM derelativization redshift; in the neutrino case, for a given HDM density, 
the derelativization redshifts are "quantized" by the discreteness of Qi,. 

By evaluating the transfer function for volatile models, we shall see that the 
points (i) and (ii) have a modest impact (see Table I in sec. 5). Consequences of 
(iii) would be important only if other cosmological parameters, like and h, 
were very well known. Altogether, therefore, the main differences induced by 
the above 3 points concern particle aspects. For example, f 's number density 
can be significantly smaller than z/'s and their masses can be much greater. In 
spite of that, their density and derelativization redshift z^er can be equal and 
LSS effects are quite similar. 

Let then still be = 7, but let us assume that the decay occurs after BBNS; 
also in this case, there are constraints related to BBNS. In fact, decays produce 
an equal amount of f 's and 7's. The present CBR temperature normalizes 7 
abundance after the thermalization of decay photons. Accordingly, a large 
V density implies a small 7 and massless u density before the decay. BENS 
would therefore take place almost at the usual temperature, but earlier in 
time. Neutron decay would be therefore allowed a smaller time to work and 
neutron abundance, at the opening of the so-called Deuterium Bottleneck, 
would be greater. From this qualitative framework quantitative limits to the 
final V abundance allowed by BENS can be derived. We shall not further detail 
this point here and will only consider a softer limit, ensuing from the obvious 
requirement that the final v density shall be smaller than 7 density, as the 
latter one has a "contribution" from "primeval" 7's. As is shown in Appendix 
A, for a model with given 0,^ ^-nd d, at redshifts z ^ Zder, the following relation 
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holds: 

Ph/p-y^^hd . (2.11) 
Henceforth, in this case, it must be 

Vthd < 1 (7 decay) (2.12) 

Therefore, according to eq.(2.8), this case never imphes primeval temporary 

matter-dominated expansion, at t < Tdy. It should be also outlined that these 
models can have quite a low massless neutrino density, as X decay increases 
7 temperature, but does not act on v background. In Appendix A we show 
that 

TJT^ = (4/ll)^/=^(l - n^d) (7 decay) (2.13) 

and a relevant variation of the massless sterile component (SMLC hereafter) 
can have a significant impact on 7^. 

On the contrary, if 0's are sterile scalars and the decay takes place well after 
BBNS, constraints are not so stringent. 0's will behave just as massless i/'s 
and these models can be characterized by a high SMLC density, which can 
also have a major impact in shaping the present LSS. 

Therefore, significant constraints on the models can arise from this effect, 
which, however, cannot be discriminated a priori. Let us recall that, in the 
absence of X decay, the ratio Pu/p-y — 0.68132(51,^/6) = Wo- X decay modifies 
it, turning gi^ into an effective value 

(16/7)(ll/4)^/^0,d (2.14) 

(see appendix A). In this case, no matter dominance occurs before Tdy if 

^hd < (1 + Wo)/2 (for sterile decay) (2-15) 

and 0's contribution to the relativistic component lowers the equivalence red- 
shift. When when d{l—flh) > i+Wp, also v's are still relativistic at equivalence. 
Accordingly, the equivalence occurs at either 

- or 10^(1-^^) .215) 

l + wo + n^d 9^ l + Wo + 2Qhd ^ ^ 

in the former and latter case, respectively. 

Mixed models involving a large volatile fraction, with late derelativization, 
arc therefore allowed only within this scenario. In what follows we shall sys- 
tematically analyze a large deal of mixed models, also with large Qh and d, 
consistent with this last picture, but we shall find that viable models with 
Q/jd > 1 are not so frequent. 
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3 Linear constraints 



Model parameters can be constrained from particle physics and/or from LSS. 
In this work we analyze a number of the latter constraints, which can be tested 
without discussing non-linear evolution. 

More in detail, we shall consider the following prescriptions: (i) We set the 
numerical constant Ax^,, in the spectrum (1.1), so that C2 is consistent with 
COBE 3(7 intervals for Qrms,ps (for different n's). 

(ii) Once the normalization is fixed at small k, we test the large k behaviour, 
first of all on the 8/i~^Mpc scale. The mass Mx, within a sphere of radius 
L = Xh'^Mpc, is 

Mx = 5.96 • W^'^nh^MQ{X/Sf . (3.1) 
Therefore, the cumulative cluster density 

CXD 

n{> M) = ^j2fK{p/M) j du[M/M{u)] exp(-'uV2) (3.2) 

Sc/ctm 

for M = 4.2/i-^ • lO^^M© is directly related to the value of (Tg. In eq. (3.2) 
we take M{u) defined so that the mass variance (evaluated with a top-hat 
window function) (Tm^u) = ^c/u; Sc values from 1.69 (Peebles, 1980) to 1.55 
were taken in figures 2 and 3. 

Let us then consider Nd — ri(> M)(100/i~^Mpc)^ for the above M value. 
Optical and X-ray observations give a value of Nd which is still not so different 
from ~ 4, as found by White at al. (1993) and Biviano et al. (1993). Both 
the observational value and its theoretical prediction are however subject to 
a number of uncertainties. For instance, observations have some problems to 
fix cluster masses (see also below). Prom the theoretical side, non-linearity 
effects and mechanism turning fluctuations into clusters cannot be said to be 
completely under control. It is also wise recalling that and feel the slope 
of the spectrum, around 8 /i^^Mpc, in a slightly different way. Henceforth, we 
shall further comment on erg values, after discussing the spectral slope. 

The above arguments tell us that, in our systematic search, it is wise to keep 
models with 1^ 10. 

(iii) Models can survive the previous test for the whole range of A-q, or for a 
part of it. Allowed A^ values were then used to evaluate the expected density 
parameters i7gas = aflbflcoii in damped Lyman a systems (for a review see 
Wolfe, 1993). Here 

Qcoii = erfc[4/V2a(M, z)], (3.3) 
where a{M, z) is the (top hat) mass variance (for mass M at redshift z) and a 
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is an efficiency parameter whicli sliould be ^ 1. More specifically, using such 
expression, we evaluated DLAS = flg^s x 10^ /a, taking z — 4.25, 5c — 1.69 
and M = 5 • IO^/i-^Mq. 

This choice of values, as well as the choice of a top-hat window function, 
is a compromise among various suggestions in the literature, supported by 
different arguments. E.g., values of Sc greater than 1.7 (Ma & Bertschinger 
1994) or as small as 1.3 (Klypin et al. 1995) were discussed, a value of M ten 
times greater was often considered and the quantitative relevance of changing 
the window function was debated (see, e.g., Borgani et al., 1996 and Pierpaoli 
et al., 1996). According to Storrie-Lombardi et al. (1995), observations give 
DLAS = 2.2 ± 0.6. Therefore, we passed models only when DLAS > 0.5. 
It must be outlined that varying this limit by a factor ~ 2 would cause only 
marginal changes for models accepted. Slight shifts of Qh or d usually cause 
significant variation of DLAS and this constraint turns out to be a fairly 
substantial one. Other tests concerning objects at large z, like early galaxies 
and QSO's, are less restrictive. 

(iv) Bulk velocities were also evaluated and compared with POTENT recon- 
structions of velocity fields. Here we shall report no detail on the procedure, 
which causes no constraint, at the 2 a level. 

On the contrary, a severe further selection arises from requiring that: (v) the 
extra-power parameter F = 7.13 • 10^^{as/a25y^^^ (here (73,25 are mass vari- 
ances on the scales i? = 8, 25 h'^Mpc) has values consistent with observations. 
Prom APM galaxies. Peacock and Dodds (1994) obtained P = 0.23 ± 0.04. In 
a more recent work, Borgani et al. (1997) give the interval 0.18-0.25 obtained 
from the Abell/ACO sample. However, owing to the high value (10/3) of the 
exponent, a shift of cxg by a small factor induces a major P displacement. For 
example, if non-linear effects are under-evaluated by just ~ 10 %, F encreases 
by ~ 37%. This effect is more likely for low P values, which can be due to 
fairly high ag values. Henceforth we kept 0.27 as top acceptable value, which is 
both the upper limit given for APM and the ~ 3 a upper limit for Abell/ACO. 
On the contrary we took 0.13 as lower limit, assuming that an underestimate 
of non-linear effects by ~ 6-8 % cannot be excluded. 

As already mentioned, the constraint on F, together with a suitable constraint 
on (Tg, implies a constraint on N^. Checks of erg, N^i and F are, therefore, 
strictly related. 

For instance, an observational value for as can be deduced from X-ray data 
on the gas temperature Tg in galaxy clusters. If clusters are substantially 
virialized and the intracluster gas is isothermal, the mass M of a cluster can 
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be obtained, once the ratio 



^ galaxy kinetic energy /mass 
gas thermal energy /mass 

is known; then, the temperature function can be converted into a mass func- 
tion, which can be fit to a Press & Schechter expression (see eq. 3.1), yielding 
the normalization of a{M) (variance as a function of the mass scale M) and 
hence erg. 

An alternative possibility, of course, amounts to deducing masses from galaxy 
velocities obtained from optical data. Both pattern imply some problem. 

The value of /3 is currently obtained from numerical models. Henry & Arnaud 
(1991), assuming /? = 1.2, estimated (jg = 0.59 ± 0.02 from a complete X- 
ray fiux-limited sample of 25 clusters they compiled. Various authors followed 
analogous patterns (see, e.g.. White et al. 1993, Viana & Liddle 1996). Eke 
et al. (1996), adding observational uncertainties and /3 error, claimed that 
(Tg = 0.50 ± 0.04. An essential issue to obtain such result is that Navarro et 
al. (1995) simulations allow to take P — 1 and an error ^6%. 

However, there seems to be a conflict between expected galaxy velocities, 
obtained assuming P — 1 ± 0.06 and optical data. In fact, the latter give a 
virial velocity dispersion in clusters ^ 800 km/s, a value consistent with 
(3Z 1-5 (Zabludoff et al. 1990, Girardi et al. 1993). 

It is possible that the critical assumption is that clusters are isothermal. In 
many clusters, cooling flows may play an important role and a fit to a cooling- 
flow cluster with a simple isothermal model may yield a mean emission- 
weighted temperature significantly reduced in respect to the virial value (see, 
e.g., Allen & Fabian 1998). However, according to Eke et al. (1998), cooling 
flows would cause a z-trend in disagreement with available data. More data 
on high- 2; clusters arc however needed to strengthen this statement. Other au- 
thors (Frcnk et al. 1990, Borgani et al. 1997) claimed that the conflict between 
predicted and observed originates from contamination of optical data by 
groups accreting onto the clusters. 

We can conclude that models with erg in an interval 0.45-0.75 should be vi- 
able. As expected, models passing previous tests keep comfortably inside such 
interval. The top Ug needed by models in flg. lb is 0.70 (for il^ = 0.26, d = 1, 
n = 1.2) and the bottom (jg is 0.50 (for Qh = 0.18, d = 16, n = 1.2). If the 
constraint on F is dropped (fig. la), just two models with ug = 0.70 (both 
yielding F < 0.10) and a model with erg = 0.48 {fth = 0.22, d = 16, n = 1.3, 
F — 0.11) are added. 

Some of the transferred spectra we obtain will be plotted, together with spec- 
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tral data obtained from Las Campanas survey (LCRS), kindly provided by 
Lin et al. (1996). It is however worth outhning soon that such comparison, 
although suggestive, is not so discriminatory. 

A detailed description of LCRS is given by Schectman et al. (1996). The survey 
encloses 3+3 "slices" 1°.5 x 80° wide, in the northern and southern emispheres. 

Data were taken using two multifiber systems. Fields 1°.5 x 1°.5 wide were 
inspected above suitable photometric limits, chosen so that there were more 
galaxies per field than available fibers. Then, target galaxies in each field were 
randomly selected and a "field sampling factor" / was memorized, to be used 
in any further statistical analysis. The average values of / are different for 
the two multifiber systems, which are able to inspect 50 and 112 objects, and 
amount to 0.58 and 0.70, respectively. The nominal photometric limits are also 
different for the two systems, amounting to 16 < m < 17.3 and 15 < m < 17.7, 
respectively. A further geometric effect is due to the impossibility to inspect 
galaxies, in a given field, if nearer than 55" . 

The actual situation for power spectrum measurement from LCRS appears to 
fall into two regimes. On scales L < 80-100 /?r^Mpc (k > 0.2 h) a fair deter- 
mination of the spectrum is obtained. In this range, LCRS results strengthen 
results from other surveys. Larger scales would be more discriminatory, but 
here errors are greater and the sample variance might cause further shifts. 

When compared with such observational spectra, a reasonable transferred 
spectrum should lie below observational errorbars, up to ~ 80-100 /i~^Mpc; 
the gap between theoretical and observational spectra is related to the amount 
of bias (constant gap means a scale-independent bias level) . For k values much 
beyond this scale, requiring a detailed fit may be excessive. Although there are 
mixed models which provide it, we see no reason to disregard models whose 
spectrum falls within ^ 3a errorbars from the reconstructed spectrum. An 
example of such models are ACDM models, with various A contents, which 
were shown also by Lin et al. (1996), but not treated in this article. 

The basic issue following this discussion, which can be suitably rephrased 
for other data sets, is that it is non-trivial to extract discriminatory criteria 
from observational spectra. Most of such criteria are expressed by quantities 
already introduced earlier in this section. Furthermore, quite in general, in 
a systematic search, predictions are to be fit with fairly wide observational 
intervals. In fact, a marginal agreement with the data, for a model met in 
a systematic analysis, can be often improved by a suitable finer tuning of 
parameters. We verified that, in several cases, this was the actual situation. 
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4 Inflationary models yielding n > 1 



This section is a quick review of results in the literature, aiming to show that 
there is a wide class of inflationary models which predict n > 1, but ^ 1.4. It 
can be therefore skipped by those who are aware of such results. 

Density perturbations arise during inflation because of quantum fluctuations 
of the inflaton fleld (f. Their amplitude and power spectrum are related to the 
Hubble parameter H during inflation and to the speed (p of the slow-rolling- 
down process, along the scalar field potential. It can be shown that the critical 
quantity is the ratio W{k) = H'^/(p, where H and (p are taken at the time when 
the scale 27r//c leaves the event horizon. The value of the spectral index can 
then be shown to be 

a(log k) 

and, if W (slowly) decreases with time, we have the standard case of n 
(slightly) below unity (it should be reminded that greater scales flow out of 
the horizon at earlier times). This decrease is due to an acceleration of the 
downhill motion of p> and an opposite behaviour occurs if (/? decreases while 
approaching a minimum. The basic reason why a potential yielding such a 
behaviour seems unappealing, is that the very last stages of inflation should 
rather see a signiflcant (/7-fleld acceleration, ending up into a regime of damped 
oscillations around the true vacuum, when reheating occurs. 

These objections can be overcame if the usual perspective is reversed: reheating 
does not arise when an initially smooth acceleration flnally grows faster and 
faster, as the slope of the potential steepens; on the contrary, reheating starts 
abruptly, thanks to a first order phase transition, perhaps to be identified with 
the break of the GUT symmetry. Before such transition and since the Planck 
time, most energy content naturally resided in potential terms, so granting a 
vacuum-dominated expansion. This picture of the early stage of the cosmic 
expansion is the so-called hybrid inflation, initially proposed by Linde (1991a). 

A toy-model realizing such scenario (Linde 1991b, 1994) is obtained from the 
potential 

V{<p, x) = (/^' - Ax')' + '2gVx' + mV (4.2) 
depending on the two scalar fields </? and x, expected to be slowly and lastly 
evolving, respectively. If the slowly evolving field is embedded in mass terms, 
the potential reads 

V{x) = M'x' + Ax' + (4.3) 

where 

A^ = / + mV' and = 2(/(^' - A/x') . (4.4) 

Eq. (4.3) shows that V has a minimum at x = 0, provided that > 0. If 
< 0, instead, the minimum is for x= ^— M2/2A, yielding /i when (p — 0. 
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Large </? values therefore require that x vanishes and then the potential 

V{ip, 0) = / + mV^ (4.5) 

allows a Planck-time inflation, as (p rolls downhill taking V from an initial 
value ~ tpf to ~ fi^. This downhill motion is expected to decelerate when the 
second term at the r.h.s. of eq. (4.5) becomes negligible in respect to fi'^, which 
essentially acts as a cosmological constant. This deceleration abruptly breaks 
down when the critical value (pc — V^pf'/g is attained, for which changes 
sign. At that point the configuration x = is unstable and the transition to 
the true vacuum configuration x reheats (or heats) the Universe. 

There are a number of constraints to the above picture, due to the require- 
ments that at least 60 e-foldings occur with (p > (p^ and that fluctuations have 
a fair amplitude. Such constraints are discussed in several papers (see, e.g., 
Copeland et al. 1994, and references therein) and cause the restriction 1.4. 
Henceforth, in this work we shall debate models with 1 < n < 1.4. 

Let us however outline that hybrid inflation is not just one of the many possible 
variations on the inflationary theme. In spite of the apparent complication of 
the above scheme, it is an intrinsically simple picture and one of the few 
patterns which can allow to recover a joint particle-astrophysical picture of 
the very early Universe, as was naturally hoped before it became clear that 
the Higgs' held of the GUT transition could not be the inflaton (see, e.g., Shafl 
& Vilenkin, 1984 and Pi, 1984). 



5 Transfer functions and spectra 

Any realistic cosmological model is expected to contain baryons (density pb — 
^bPcr), photons (present temperature To = 2.733^ K, density p^), CDM (den- 
sity Pc — ^cPcr) and a SMLC. The basic SMLC are the usual 3 massive i/'s. 

Mixed models include also a HDM component (density ph = ^hPcr) and the 
SMLC is however modifled. If HDM particles are massive i/'s, SMLC density 
is lower, or even vanishing when all z/'s are massive (of course, sterile massless 
particles can be however added ad hoc). In the volatile case, we face the 
opposite situation, as the SMLC must contain at least a scalar, in top of 
usual massless i/'s. In all cases, we can parametrize SMLC by using w — 
0.68132(^^,e///6), according to eq. (2.14). 

The algorithm used to evaluate the transfer functions of a wide range of mixed 
models is discussed in Appendix B and, besides of allowing any value of fib, 
flh, w and other standard parameters, allows to flx arbitrarily the expression 
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of $(j9) [sec cqs. (2.1) and (2.7)] and the dependence on p of the amphtude of 
HDM fluctuations. 

As aheady mentioned, however, in this work Q = 1,A = 0, i? = 50 km s~^Mpc~ 
and Qfe = 0.1 (some plots for thermal models, with different values of are 
however shown). Exploring other values of the former 2 parameters would 
amount to extending the physical range of models. The expected behaviour 
for different values of the latter two parameters, instead, can be qualitatively 
inferred from the results shown below. A test of the effect of varying them 
should be focused onto a portion of the parametric space, selected on the 
basis of the results of this work and other physical criteria. 

Within the above restrictions, we tested models at regular logarithmic interval 
for Zder-i taking log2 d = — 1,...,5 (7 values) and at regular interval for hot- 
dark-mattcr density, taking VLh = 0.10, 0.12, 0.44, and 0.45 (19 values). 
A non-systematic sampling of this parameter space allowed us to exclude some 
models a priori, reducing the total exploration to 120 models. 

There is a significant overlap between these values of Qh and d and those 
allowed in mixed models with massive z/'s, as can be seen also in fig. 1 here- 
below. One of our basic aims was then to obtain a quantitative estimate of 
the impact of the different momentum distribution on transfer functions. To 
this aim we re-evaluated T, for a few mixed models with HDM made either 
by massive neutrinos or volatiles, but with the same values of Q^, d and gu,eff, 
according to eqs. (2.3), (2.4), (2.14). The results of this comparison for two 
model sets are shown in Table I. 

The detection of the (modest) amount of such shifts is one of the results of 
this work. However, this does not decrease the relevance of volatile models, 
as they allow parameter choices which, otherwise, would be impossible with 
the standard 3 //-flavours. Furthermore, these models require and allow a high 
level of SMLC. 

The main output of this work is however represented by fig. l(a,b), where 
we show which parts of the parameter space are compatible with linear con- 
straints, for spectral indices n = 1, 1.1, 1.4. 

Fig. lb takes into account all the constraints mentioned in sec. 3. Models inside 
the curves fulfill the requirements on COBE quadrupole value (3(7 intervals), 
galaxy clusters, high-z objects and F. On the contrary, fig. la cxchidcs the 
requirement on F. Models forbidden by fig. lb but allowed by fig. la should 
be accurately tested for F in non-linear simulations. 

For most models considered, the spectrum 

p{k)^Ak''r\k), (5.1) 
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Table 1 

Ratio between tranfer functions in volatile and thermal models (T (vol)/ T (ther)) 
at various comoving scales for different components. The upper part refers to models 
with = 0.3 and 3 massive neutrinos, the lower part to models with 0,^ = 0.1 and 
1 massive neutrino. In both cases CI = 1, h = 0.5 and flf, = 0.1. 



L/Mpc 




cold 


hot 


baryons 




total 


0.2000E+05 




1 


1 


1 




1 


0.5120E+03 




1 


1 


1 




1 


0.2048E+03 


0. 


,9995E+00 


0.9996E+00 


0.9995E+00 


0. 


,9995E+00 


0.8192E+02 


0. 


,1001E+01 


0.1002E+01 


O.lOOlE+01 


0. 


,1001E+01 


0.3277E+02 


0. 


,1003E+01 


0.1004E+01 


0.1003E+01 


0. 


,1003E+01 


0.1311E+02 


0. 


,1003E+01 


0.1004E+01 


0.1003E+01 


0, 


,1003E+01 


0.5243E+01 


0. 


,1005E+01 


0.1003E+01 


0.1005E+01 


0. 


,1004E+01 


0.2097E+01 


0. 


,1007E+01 


0.1012E+01 


0.1007E+01 


0. 


,1009E+01 


0.8389E+00 


0. 


,1006E+01 


0.1036E+01 


0.1006E+01 


0. 


,1015E+01 


0.2000E+05 




1 


1 


1 




1 


0.5120E+03 




1 


1 


1 




1 


0.2048E+03 


0. 


,9995E+00 


0.9996E+00 


0.9995E+00 


0, 


,9995E+00 


0.8192E+02 


0. 


lOOlE+01 


O.lOOlE+01 


O.lOOlE+01 


0. 


,1001E+01 


0.3277E+02 


0. 


,1003E+01 


0.1002E+01 


0.1003E+01 


0. 


,1003E+01 


0.1311E+02 


0. 


,1003E+01 


0.1002E+01 


0.1003E+01 


0. 


,1003E+01 


0.5243E+01 


0. 


,1005E+01 


O.lOOlE+01 


0.1005E+01 


0. 


,1005E+01 


0.2097E+01 


0. 


,1007E+01 


0.1008E+01 


0.1007E+01 


0. 


,1007E+01 


0.8389E+00 


0. 


,1006E+01 


0.1026E+01 


0.1006E+01 


0. 


,1008E+01 



where the constant 

A = ^-^x-o^'A^, (5.2) 

starts - at low A;'s - with values similar to standard CDM. If n > 1, it soon 
abandons CDM behaviour, raising in a slightly steeper way. Both for n > 1 
and n = 1, however, its bending at maximum is sharper. At the r.h.s. of 
the maximum it returns below CDM. At large /c's (~ lMpc~^), its decrease, 
sometimes, is (slightly) less steep than CDM. 

Some examples of such behaviour are given in fig. 2 for volatile models and in 
fig. 3 for thermal models; all cases we show have n > 1. Models are compared 
with CDM (dotted curve) and with the spectrum reconstructed by Lin et al 
(1996) from LCRS data (3(7 errorbars). Models are ordered with increasing 
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Qh- For volatile models, various values of d arc considered, including high 
values; in particular, in fig. 2c we show a spectrum obtained with Zder — 625. 

Massive neutrino models obviously have w < Wg. Extra SMLC can be however 
added ad hoc and the comparison of Table I was made after implementing w 

to the value of volatile models. The boundaries shown in fig. 1 are also valid 
when all SMLC expected in a volatile model is included. Its presence causes a 
later transition from radiation dominated to matter dominated expansion. As 
a consequence, the peak of the transferred spectrum is shifted to smaller /c's 
and, in general, this favours the agreement of models with linear constraints. 

While the boundaries marked in fig. 1 therefore apply to neutrino models with 
extra SMLC, the spectra shown in fig. 3(a-b) are for physical models without 
extra SMLC. The spectrum in fig. 3c, instead, is obtained adding an amount 
of SMLC corresponding to 3 massless i/'s. 

The inverse transfer functions of all models considered, were also parametrized 
with 4-th degree polynomials in y/k, both ior z — and z — 4.25. The table 
of coefficients will be made available upon request. 



6 Discussion 

After evaluating the transfer function of a set of more than 120 mixed models, 
with Q — 1 and h — 0.5, we obtained the expected values of several observable 
quantities, which can be estimated using the linear theory. In particular we 
worked out the expected cluster number density, the abundance of Damped 
Lyman a clouds, the mass variances at 25 and 8 /i~^Mpc. From them we also 
estimated the extra power parameter F. 

The HDM of the models we considered was cither thermal (massive i/'s) or 
volatile (arising from heavier particle decay). We discussed the parameter 
space for the latter case and compared the transfer functions for thermal and 
volatile HDM, for values of flh and d which allow both models. Differences 
between thermal and volatile models are then really significant only as far as 
the expected SMLC is concerned. 

Results on parameter constraints are summarized in figures 1 (a,b). As was 
already known, a number of such models with 0,^ up to 0.30 pass the above 
tests. If n— 1, volatile models allow little extra freedom, namely for high z^er- 
The situation is already different for n ~ 1.1. Here models with Zder — 600 
are allowed for Qh up to 0.14 and greater Qh are allowed for values of z^er still 
much lower that those allowed by neutrino models. The range of flh values 
allowed with z^er — 600 extends upwards as n increases and overcomes 0.20 for 
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n = 1.4. The greatest value found for Q/j is 0.30 for Zder = 10'* with n = 1.2. 
Thermal models with 2 or 3 massive neutrinos and a suitably added SMLC 
arrive to ilh ~ 0.28, for n = 1.2 and n = 1.3, respectively. 

These values may not seem too large, in respect to Qh values currently used 
in the literature (with thermal models). It must be outlined, however, that 
our acceptance criteria are more stringent than usual. With such criteria, no 
thermal model with flh > 0.16 is accepted for n < 1. 

Let us also draw the attention on the very low values of z^er that volatile 
models with large n allow. As an example of low Zder, in fig- 2c we show the 
transferred spectrum of a model with 22 % of HDM and Zder = 625, which, 
for a fairly high value of COBE quadrupole, has excellent fits with all linear 
constraints. HDM particles of this and similar models have a mean square 
velocity ~ 500km/s today, still away from potential wells. Such speed is likely 
to guarantee them not to cluster with baryon or CDM on any scale. Dynamical 
mass estimates, in a world containing such component, might lead to observe 
~ 75-80 % of critical density. 

It should be also mentioned that, if the constraint on F is dropped, greater 
QhS seem allowed (see fig. la). The top value we found is Qh = 0.45 with 
Zder — 5000. With 3 or 2 massive neutrinos the highest ilh obtainable are 0.40 
and 0.38, respectively. As already outlined, in fig. 1 neutrino model curves are 
overlapped to volatile model parameter space, but the allowed regions apply 
to them only after adding suitable SMLC. 

Liddle et al. (1996), in a paper focused on n < 1 models, mentioned that a 
model with Qh — 0.35, n — 1.2, h — 0.4 agreed with the finear constraints they 
imposed. Using our linear constraints, is to be lowered to ~ 0.30. Volatile 
models do not need to lower h so far, to agree with data, as they naturally 
have SMLC, which induces fairly similar effects. Also thermal models can 
agree with data keeping to h — 0.5, rising n, with or without extra SMLC. 
Fig. 3 illustrates some spectra for such models. 

Models with high HDM content and high average kinetic energy, could also 
ease the problem of the apparent baryon excess in several galaxy clusters. It 
has been known for several years (see, e.g.. White & Frenk 1991, Boehringer 
et al. 1992, Briel et al. 1992, David et al. 1993, White et al. 1993) that the 
baryon mass fraction in galaxy clusters is high, in respect to expectations in 
flat, pure CDM models. According to White et al. (1993), even neglecting 
baryons in galaxies, the observed baryon fraction yields 

Qfe/i^/^ = 0.05 ± 0.02 (6.1) 

for fiat, pure CDM models. BENS limits of Walker et al (1991) required then 
that flbh'^ = 0.0125 ± 0.0025. 
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Much work was then performed to analyze individual clusters and lists of 
baryon contents in cluster samples were presented by several authors (see, 
e.g., White & Fabian, 1995). The situation is summarized in a recent paper 
by Evrard (1997) who obtains the constraint 

n^/i3/2 = 0.060 ± 0.003 (6.2) 

This is to be taken with recent BBNS hmits on Qi,- E.g., Copi et al (1995) 
give 0.007 < Qft/i^ < 0.024 and therefore ~ 0.1, with h = 0.5 is marginally 
allowed. For h = 0.5, as wc assumed here, eq. (6.2) gives a central value 
Qb = 0.17 and a 3cr minimum Qb — 0.14. 

If HDM does not cluster with CDM and baryons, a measure of Qb would 
however yield 0,b,app = ^b/i^^^h)- Within mixed models, therefore, agreement 
with data requires > 0.41 or > 0.29, respectively. We have shown in this 
work that the latter value is coherent with LSS constraints, if 1.2 < n < 1.4; as 
n increases, compatible models have Zder values decreasing, down to ~ 2500. 

Limits become more permissive if masses obtained from lensing arc replaced to 
those obtained from X-ray data (see, e.g., Allen & Fabian 1998). High baryon 
density, however, seems a widespread feature. Perhaps the limits shown in 
eq. (6.2) are to be reset, taking greater cluster masses, but, if = 1, there 
seems to be a definite evidence of a component which tends not to cluster with 
CDM. The debate is related to the discussion on as hmits in sec. 3. 

In the literature, mixed models with n < 1 were often considered and allow to 
predict acceptable values of a number of observable quantities, if low rates for 
the HDM component arc taken. Considering n > 1 is at least as legitimate as 
taking n < 1 and leads to a range of mixed models allowing fair predictions 
on the same quantities. We showed that models with n > 1 require and allow 
higher HDM contents. They could therefore improve our understanding of 
why fto measures seem to give increasing values when greater scales L are 
considered. This was one of the basic motivations to introduce mixed models, 
more than a decade ago. 
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FIGURE CAPTIONS 



Fig. 1 - The curves enclose areas where models predict LSS linear features 

consistent with observations. Each curve refers to the value of n written aside. 
Point-dashed lines refer to models with HDM made of massive z/'s. See the 
text for a more detailed comparison of the LSS they predict with volatile 
models. In fig. la the constraints taken into account are: cluster abundance, 
(Tg value, gas in damped Ly-a systems. In fig. lb the value of the extra-power 
parameter P is also constrained. 

Fig. 2 - Mixed model spectra with volatile hot component (solid line), com- 
pared with standard CDM (dotted line). The LCRS reconstructed spectrum is 
also reported (Scr errorbars). Models a, b, c are consistent with all constraints. 
Model d is outside the contours plotted in fig. lb and has a small P. Models 
like it would agree with data for a still greater n. 

Fig. 3 - Mixed model spectra, with thermal hot component (solid line) com- 
pared with standard CDM (dotted line) and LCRS reconstructed spectrum, 
as in fig. 2. Models a.b agree with all constraints. Model c is in marginal 
disagreement, in spite of adding extra SMLC amounting to 3 massless z/'s. 
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Appendix A 



In the decay X ^ v + (p, both v and (p acquire momentum (and energy) 
Pin = If a decay occurs at tp, at any later t, the momentum is redshifted 

into p — Pina{tp)/a{t); in turn this means that, at any t, a t'-particle of 
momentum p was born when the scale factor was 

a{tp) = a(t) Iplmx ■ [oX) 

The map of p on a(ip) yields that 2a{t)/mx — chitpj/p and 

dp/d[a{tp)] = mx/Mt) = PKtp) ■ (a2) 

This allows to give the momentum distribution of f-particles 

^ J_^^ = ^MM^ exp {-^\ (a3) 

Nv dp pTdy V ^dyj 

when the decay is complete. In fact, using (a2), dNy/dp —[dNy/dtp] [dtp/ da{tp)] 
[da{tp)/dp] =[N^/Tdy\[l/a(tp)][a(tp)/p\. 

The expression (a3) is greatly simplified is the expansion is radiation domi- 
nated during the decay. In such a case a{tp)/a{tp) = 2tp, while, using eq. (al), 

Q = tp/Tdy = [a(tp)/ady]'^ = (p/py (a4) 

with 

p^ {'mx/2)[ady/a{t)] . (a5) 

Accordingly, 

$,(p) = 2(Q/p)exp(-g) . (a6) 
Henceforth, at large p, the distribution decays oc exp(— p^). Large p come from 
late decays, when the decay process is however radiation dominated; therefore 
the latter property is always true. 
Owing to eq. (a4) 

dQ/Q^2dp/p (a7) 
and, therefore, can be easily verified to be normalized to unity. In a similar 
way 

oo oo 

(p) = J dpp^yip) ^"^P j dxx^ exp{-x^) ^p^j^ (a8) 



This allows to replace, in the definition (a5) of p and, henceforth, in the 
distribution parameters referring to the decay process (rax and Tdy) with 
parameters explicitly related to a cosmological model: the derelativization 
time/redshift /scale factor (subscript der) and the v-particle mass m. By defi- 
nition, let the derelativization occur when (p) = m or, according to eq. (a8), 
when p— 2ml \pK. Then, at any time: 

p= {2m/^/Ti){ader/a) . (a5') 
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Henceforth, the distribution peak and cut-off move towards smaller energies 
as the Universe expands, with pa = const.. 

In Appendix B, we discuss the time evolution of a model. The basic refer- 
ence redshift will be Zeq = 10"^. If ^ = To/2.7333 K= 1, the present radiation 

density Pr,o = Pcr/Zeq. 

In respect to Zeq let us define: d — Zeq/ z^eri U = Oj^q/a = (1 + z)6* /IQ'^A.h. It 
will be: 

(p) — myd and p — myd2/^/^:. (a9) 
In the case of HDM made of massive i/'s, the momentum distribution reads: 

= [2/?,a'i)]\p/TA'W'^' + 1)]-' (^6) 

and, accordingly, 

(p) = [77rVl80C(3)]T, (n8) 

Eq. (n6) shows that, for large p, such distribution decays oc exp(— p). 

Furthermore, z/'s become non-relativistic when Ty = ml80C(3)/77r^ and, 
while the former equation (a9) still holds, it is 

Ty = myd 180C(3) /77r^ . (n9) 

Let us now show that 

Ph/ P-y — ^hd at high z. (alO) 

At such z's, in fact, (e) = (p) can be expressed through eq. (a9). In turn, the 
V number density can be obtained rescaling its value at z — 0. Henceforth 



, MhPcr faoY o ^ 
Ph^myduh^myd — = Q^d 

m \ a J 



Per [ 

ao 



ao\ 
a J 



4 



(all) 



and, owing to the definition of a^q, this shows (alO). 

Let us now consider the constraints to a volatile model when = 7. The 
photon density before decay {pj^eariy) has safe lower limits related to BBNS, 
that we shall not discuss here; a fresh 7 component, whose density coincides 
with pfi, accretes onto it because of the decay; hence, at any time, p^ > ph 
and, owing to eq. (alO), 

D,hd <^ 1 for 7 decay. (al2) 

Under such condition, it is almost granted that no matter dominance occurs, 
during the decay process. The actual constraint for this not to happen is that 

flhd < [1 + 2/(1 + ^^/8.8)]~^ for 7 decay. (al3) 

In fact, comparing densities before and after r^y, we have that the top value 
of px ^ 2ph and p^.eariy - P7 " Ph] furthermore p^ = WoPj,eariy with Wo = 
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(4/11)^/^7(71^,0/16 ^ gu,o/8.8 {g,y^o are the spin states of masslcss z/'s). The 
requirement that the top value of px < P'y,eariy + Pu, yields 2ph < (l + u'o)(p^ — 
Ph) and hence, owing to eq. (alO), flhd -C [l + 2/(l + Wo)]~^. E.g., for g^^o = 6, 
eq. (al3) yields Qhd < 0.46. 

Let us now consider the case of massless sterile scalar 0's. Arguments similar 
to those leading to eq. (al3) assure no matter dominance before Tdy if 

flhd < (1 + Wo)/ 2 (for sterile decay). (^14) 

In fact, in this case, the cosmological effects of 0's add to those of massless 
i/'s and we shall therefore distinguish the real v density pi,^o from the effective 
v density 

Pu = Pu,o + P,f> = wp^ {alb) 
where w = Wo + flhd. Once again, the top value of px is — 2ph and, in order 
that, at all times, px < P'y{l + Wo), owing to eq. (alO), eq. (al4) must be 
fulfilled. The increase of the radiative sterile component can be also expressed 
renormalizing g^^o into 

9.,eff = g.,o + 8.80640 Qhd (al6) 

which coincides with eq. (2.14). 

Let us then discuss where the actual equivalence redshift Zeq is set. If v- 
particles do not contribute to the relativistic component, 

Zeq = Per/ (1 + w)p^ = Zgg/ (1 + w)6''' (al7) 

If Zder < Zeq, as is defined here, i.e. ii d{l — flh) > 1 + Wo, v^s are to be added 
to the relativistic component and 

Zeq = Zeq(l - flh) / {1 + + 20^0?)^^ . (al8) 



Appendix B 

The system of equations to be integrated in order to work out fiuctuation 
evolution has been discussed in many papers. E.g., an excellent and updated 
discussion is given by Ma & Bertschinger (1996). Here we shall report only 
the changes needed to treat HDM, if it is made of w-particles. 

Let the phase-space distribution be / = /o(l + e). For massive i/'s 

/„ = (27r)-3y,(e" + 1)-^ with x^p/T, (nl) 
and, therefore, the expression 

{x/fo)idfo/dx) = -x/{l + e-^) (n2) 
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will enter the equation (b3) here below. For volatiles, instead, 

fo = {2tx)-^ N,{ao/apfQe-^ (61) 

and 

{ap/fo){dfo/dap) = -(1 + 2Q) (62) 
will be used in the equation: 

where k is the wave-number of the perturbation, H accounts for the grav- 
itational interaction and /i yields the angle dependence of e. It should also 
be outlined that, while for massive i/'s, p/E — [1 + l/a;^|/^d^]~^/^, with — 
c^l80C(3)/77^^ here 

1/2 

(64) 



1 

1 + 



P_ ^ 

E ['^ Qy'dl 
with = d2/ y^. 

The different phase-space distribution has further effects on the gravita- 
tional field equations, that we treated in the synchronous gauge. Here we 
consider the stress-energy tensor Tjj for massive i^'s and its perturbation 

A(2roo -T)^ j d^pE{l - m'/2E^)fo<To, (n4) 

where ao is the 0-th order term in a spherical harmonic expansion of e. As 
E{1 - w?/2E'^) = m{l/2 + x^y'^dl)/{l + x^y^dlf/^ it turns out that 

A(2roo - r) = ^pAz^y^Aydr^) (^5) 

with 



oo 



e3,.(s) = 2 I dxe -n{x)-^^^^^ (n5') 







where 7r{x) — x^/{l + e For volatiles we have an analogous expression: 



A(2roo - r) = -p,Az'y'e^,,{yd,) (56) 



with 



oo 



e3,(.) = 2 / cige-«^i/^±|^ (65') 







In the gravitational field equations also the quantity 64 = /q°° dxe^^ XT:{x)ai{x) 
is used, for massive i^'s (a"i is the first order term in the spherical harmonic 
expansion of e). For volatiles it ought to be replaced by 



e4,„ = j dQe-^^ai{Q) . (66) 
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Our code integrates the set of equations for the fluctuation evolution in 
the synchronous gauge. We sped it up by truncating the hierarchical set of 
equations of free-streaming particles (photons after recombination, HDM and 
massless i/'s), according to the scheme suggested by Ma & Bertschinger (1996), 
slightly modifying it because we don't use conformal time. Results with trun- 
cation were compared with results with a free number of harmonics and we 
found that a truncation at the 8th harmonics is mostly sufficient. Results 
presented here are however obtained with truncation at the 24th harmonics. 
Models were evolved down to z — 0. We evaluated the transfer function for 
24 scales ranging from L = 20000 Mpc to L = 0.53054 Mpc. Results from 
the public code CMBFAST, yielding the transfer function for mixed models 
(only with thermal hot component), were compared with the results of our 
algorithm, for a number of models and for scales fixed by CMBFAST. 

Relative discrepancies between code outputs vary from less than 10""^, on 
scales above ~ 100 Mpc, to shghtly above 3 % in the worst cases, which occur 
for the hot component, when 0,^ > 0.3 and for scales below 2-3 Mpc. In 
such cases CMBFAST tends to give a transfer function slightly smaller than 
our code. The typical discrepancy, however, keeps safely below 0.5 % and, 
therefore, we fully con&m the validity of CMBFAST. 

Altogether the main source of error comes from the analytic fit to the 
transfer functions, although the top discrepancy at a single point keeps < 1 % 
in the worst cases and its typical value is < 0.1 %. 
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logg(d) 




logg(d) 



CDM: Q=15.7A6K,n=l.N^,=43,ng=23,r=0.42 

MDM: Q=14.2MK,n=l.l.nh=0.16,Zj„=5000 
N^,=8- 13,DLAS= 1 .2.r=0.20.a8=0.67 
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CDM: Q= 15.7MK.n= 1 .N^,=43.ng=23.r=0.42.nb=0. 1 

MDM: Q=12.3MK.n=l.l.nh=0.16.g^=6.nb=0.09 
N^,=7- 12.DLAS=3.7.r=0.23.a8=0.65 
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log(k Mpc/h) 



CDM: Q= 15.7A6K.n= l,N^,=43.ng=23.r=0.42 

MDM: Q= 16.0MK,n= 1 .2.0^=0. 18,Zdep= 1250 

N^j=2.3-4.4,DLAS=0.7.r=0.15.a8=0.56 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2b - 




log(k Mpc/h) 



CDM: Q= 15.7MK.n= 1 .N^,=43.ng=23.r=0.42.nb=0. 1 

MDM: Q= 12.6MK.n= 1.2.nh=0.27.g^=6.nb=0.09 
N^,=6-10.DLAS=0.8.r=0.16.a8=0.65 




-2 -1 

log(k Mpc/h) 



CDM: Q=15.7MK,n=l,N^,=43.ng=23,r=0.42 

MDM: Q=12.6MK.n=1.3.nh=0.33.g^=4 (+SMLC) 
N^,= 17-23.DLAS=0.5.r=0.14,CT8=0.80 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

3c - 




log(k Mpc/h) 



CDM: Q= 15.7A6K.n= l,N^,=43.ng=23.r=0.42 

MDM: Q=11.3MK,n=1.4.nh=0.30,Zdep=2500 

N^,=3.2-5.7.DLAS=0.7.r=0. 1 1 .a8=0.59 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2d - 




log(k Mpc/h) 



CDM: Q= 15.7A6K.n= l,N^,=43.ng=23.r=0.42 

MDM: Q=16.lMK,n=1.4.nh=0.22,Zdep=625 

N^,= 1.3-3.6.DLAS=2.8.r=0.13.a8=0.54 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2c - 




log(k Mpc/h) 



